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We construct an effective Hamiltonian for the motion of T^g highly correlated states in Na^CoC^. 
We solve exactly a multiband model in a C0O6 cluster with electronic occupation corresponding 
to a nominal Co valence of either +3 or +4. Using the ensuing ground states, we calculate the 
effective O mediated hopping t = 0.10 eV between many-body T2 g states, and estimate the direct 
hopping t' ~ 0.04 eV. The trigonal splitting 3D = 0.315 eV is taken from recent quantum chemistry 
calculations. The resulting effective Hamiltonian is solved using a generalized slave-boson mean-field 
approximation. The results show a significant band renormalization and a Fermi surface topology 
that agrees with experiment, in contrast to predictions using the local-density approximation. 

PACS numbers: 71.27.+a, 71.18.+y, 74.70.-b, 74.25.Jb 



The doped layered hexagonal cobaltates Na^OoC^ 
have attracted considerable interest in the last few years. 
The exceptionally high thermopower and at the same 
time low thermal conductivity and resistivity for 0.5 < 
x < 0.9 Q, 0, makes the system attractive because of 
its potential technological applications. Wider attention 
has been triggered by the discovery of superconductiv- 
ity in hydrated Nao.asCoC^— 1>3 H2O 3|. Because of the 
hexagonal structure of the system which frustrates com- 
peting magnetic interactions, a possible natural explana- 
tion of the superconductivity seemed to be a resonance- 
valence-bond state in an effective one-band model |4j. 
However, based on geometrical and physical arguments, 
Koshibae and Maekawa pointed out that it is more real- 
istic to start from a three-band model of degenerate Co 
3d t2g orbitals with an indirect hopping over intermedi- 
ate O orbitals This leads to four independent in- 
terpenetrating Kagome sublattices. On the other hand, 
local-density approximation (LDA) calculations predict 
a Fermi surface (FS) with a large portion around the 
r point of the Brillouin zone, and six pockets near the 
K points 0, II] ■ However, angle-resolved photoemission 
(ARPES) revealed onl y a single hole like FS around F 
and no hole pockets [9j, llOi lll| . Recent experiments with 
X-ray absorption spectroscopy (XAS) and their 
interpretation |12lll3lll4| show a strong Co-0 covalency. 
In particular, the amount of the 3d n configuration in a 
C0O6 cluster is smaller than the total amount of 3d n+1 L 
states (L denotes a ligand O 2p hole) for both, nomi- 
nal Co 3+ (n = 6) and Co 4+ (n = 5). Ab initio quantum 
chemical calculations also obtain a large Co-0 hybridiza- 



tion 0. 

It has been suggested that the discrepancy between 
the LDA calculations and the Fermi surface measured by 
ARPES can be explained as an effect of the correlations 
[lfij |. which are absent in LDA and LDA+U approaches, 
even in the most sophisticated treatments |l7| . How- 
ever, the effective Hamiltonian used is unjustified and 
differs from that we have derived (see below). Moreover, 
the dispersion relation along T-K differs from the exper- 
imentally observed one, and the results contradict those 
obtained using dynamical mean field theory, which ac- 
tually obtain an enhancement of the hole pockets [lij] . 
Another recently proposed explanation is a localization 
of the hole pockets by disorder 0] . However this local- 
ization is not expected to alter the small volume enclosed 
by the main portion of the LDA Fermi surface around _T, 
while according to the more recent ARPES experiments, 
this portion contains all holes in a way consistent with 
Luttinger theorem 

This discussion shows the need of further studies of the 
electronic structure of Naa,Co02, and the starting point 
should take into account the strong Co-0 covalency. This 
is the purpose of this Letter. We start from the exact so- 
lution of the appropriate multiband Hamiltonian for Co 
3d and O 2p holes in the basic CoOe octahedron, with the 
hole content corresponding to x = 1 (nominal Co 3+ ) and 
with one additional hole (nominal Co 4+ as for x = 0). 
The parameters of the multiband model are fixed by a 
previous fitting of the polarization dependent XAS spec- 
tra [ljj . The resulting ground states are mapped onto the 
corresponding states of an effective three-band model for 
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(fictitious) 3d t^g orbitals on Co sites. The eigenstates of 
the cluster are used to calculate the hopping in the effec- 
tive model. The procedure has some similarities to the 
derivation of an effective t—J model in the cuprates using 
non-orthogonal Zhang- Rice singlets |2jJ, 0|- The trig- 
onal splitting is taken from recent difference-dedicated 
configuration interaction (DDCI) calculations 0] in a 
distorted CoOg cluster. The DDCI are quantum chem- 
istry calculations, particularly suited to calculate excita- 
tion energies in strongly correlated systems Ei ] . The 
effective model is solved in a slave-boson mean-field ap- 
proximation [53, 13 ■ The results reproduce the essential 
results of ARPES, in particular the shape of the FS and 
the quasiparticle dispersion near the Fermi energy ep. 

It might seem that the procedure mentioned above, 
which eliminates the O 2p states from the low-energy 
physics, is not valid when the Co-0 hopping is larger or 
of the order of the Co and O on-site energy difference A. 
However, calculations in the cuprates [25L E?3 | and other 
systems |27j using the cell perturbation method showed 
that this is not the case. In fact as the hopping increases, 
the d states are replaced by hybrid states of the same 
symmetry which remain well separated from the other 
states. Using the mapping procedure it is even possible 
to calculate the O spectral density at low energies with 
the effective model which does not contain explicitly O 
states 1 2 Cil . 

The Hamiltonian for the C0O6 cluster contains Co 3d 
and O 2p holes and includes the splitting between 3d e g 
and t2 g orbitals, Co-0 and 0-0 hoppings and all inter- 
actions between 3d holes. Spin-orbit interaction is ne- 
glected. Details are given in Ref. 0|. We assume oc- 
tahedral Oh symmetry. This assumption greatly reduces 
the size of the matrices to be diagonalized and has a very 
small effect on the resulting many-body eigenstates |l ll . 
The ground state of the cluster with four holes (nomi- 
nal Co 3+ ) \Ai g ) is a singlet with A\ g symmetry. We use 
capital letters to denote the irreducible representations 
of many-body states. For the parameters that fit the 
polarization dependent XAS spectra, it has been found 
that \Ai g ) consists mainly of 30% of the state with four 
3d e g holes and 47% of states with three 3d e g holes 
and one 2p hole, in a linear combination of e g symme- 
try [l4|. For one additional hole, the ground state of the 
cluster is a spin doublet with Ti g symmetry (six-fold de- 
generate). These states will be denoted as I7, er), where 
7 = xy,yz or zx and a =| or j. As a first approxima- 
tion, they can be constructed adding a 3d ti g hole to 

\Mg). 

In this work, to represent the movement of holes in 
the system between different CoOg clusters, we map the 
state |z,^4i g ) where i is a site index, onto the vacuum 
(absence of holes) at site i, and the states 1 1,7,0") onto 
the corresponding states /i| 7(T |0) of the effective three- 
band model, where /ij creates a hole at site i, orbital 7 



and spin a. The most important effective hopping is the 
one mediated by O orbitals. Geometrical considerations 
show that a hole in the state \i, 7, a) can hop only to some 
nearest Co atoms (building the state \j, 7', a)) through an 
intermediate O 2p orbital. Fig. ^ illustrates this indirect 
hopping, which can be well approximated by |2l| 



-t = (j, zx,t\(i, A lg \H hop \i,yz,1)\j, A lg ) 



(1) 



Hhop is the sum of Co-0 and 0-0 hopping terms of both 
clusters. For the parameters that fit the XAS spectra 
we obtain t = 0.10 eV. This agrees with previous es- 
timations Q . Considering all possible indirect hoppings 
leads to the picture of four interpenetrating Kagome sub- 
lattices for the degenerate I7, a) orbitals (see Refs. 
and Fig. [J). 




FIG. 1: (Color online) Orbitals involved in the effective hop- 
ping of a hole with symmetry yz to a neighboring state with 
symmetry zx, through a p z orbital of an intermediate O site. 

While t is the more important hopping process [|J , the 
direct hopping t\ between d electrons might be important 
for a realistic description of the band structure 0, 0] . 
Unfortunately, since this hopping is not included in the 
calculation of the C0O6 cluster used to fit the XAS spec- 
tra, it cannot be extracted using this experimental in- 
formation. Taking the value t\ — —44.6 meV obtained 
from a fit of the band structure in Ref. Jr|, we obtain 
t 1 = 36 meV > for the effective hopping between highly 
correlated hole states 

*' = (j^y^KhAigl-tidl^di^l^xy^^j^ig) (2) 

Here i and j refer to two nearest-neighbor clusters in the 
(1,1,0) direction. 

In presence of a trigonal distortion, the point sym- 
metry is reduced to D 3 d and the T 2g states are split 
into A' lg and E' g states. The prime denotes irreducible 
representations of D^d- This splitting, which is an es- 
sential parameter of the effective model, is beyond the 
reach of our cluster calculation in Oh symmetry. For- 
tunately, recent accurate DDCI calculations give a split- 
ting 3D — 0.315 eV between the many-body ground state 
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(\i,xy, a) + \i,yz, o) + \i, zx, cr))/v3 of A' lg symmetry and 
the excited E' g states. The large value of D compared to 
estimations using a point charge model || is due to the 
effects of the Co-0 hybridization. 

For D ^S> t, the one-band model is justified, while for 
D -C t, the picture based on four Kagome sublattices is 
a better starting point (Refs. 0,0 take D = 0). As ex- 
plained above, the cluster calculations give D ~ t, which 
corresponds to an intermediate situation. This is also 
supported by an analysis of the XAS intensities under 
different polarizations of the incident light [l4|. I n what 
follows we take D = t, t' = 0.4 t with t ~ 0.1 eV. Note 
that the indirect hopping of holes through O orbitals —t 
and the direct hopping t' have different signs, but we 
define them in such a way that t > 0, t' > 0. 

The energies t, D and t' are smaller by at least one 
order of magnitude in comparison with all other energy 
scales in the problem, like Co-0 charge-transfer energy 
A ~ 4 eV and Coulomb repulsion at the same d orbital 
U ~ 4.5 eV [HE1E3. Therefore, we assume that the 
energy cost to bring two additional holes in the same 
cluster is much larger than the effective hopping. This 
implies that double occupancy is forbidden in the effec- 
tive low-energy model. 

This model can thus be written in the form 

Heff = Hi + H int 

i,j 7.7'," (3) 

with U — ► oo to avoid double occupancy. Fourier trans- 
forming the hole operators, Hi takes the explicit form 

Hl=J2 E ( e 77'W+4 7 '( fe )+^77')47>fe7'- 
k 7,7',ct 

(0 COS $3 COS 02 N 
COS #3 COS 9i 
COS 02 COS 0i 

/cos 6»i 
e'^,(k) = 2t'\ cos0 2 

V cos 3/ 



D 77 , = -D I 1 1 




(4) 



and 3 = **2 -i- 



With 01 = fcjS, 02 = 2^ 2^ auu 173 — 1 

where here k x refers to the direction of two nearest- 
neighbor Co atoms ((—1,1,0) in Fig-QJ. We believe that 
H e f f is the starting point for the description of the low- 
energy physics of the cobaltates. While it is similar to 
previous proposals 0, 0] , the most important parame- 
ters D and t are derived from cluster calculations that 
properly take into account the important effects of Co- 
O covalency and all correlations inside the Co 3d shell. 



Due to the large Co-0 hybridization, the effect of the 
distortion of the CoOg octahedra is stronger than in pre- 
vious multiband approaches, but not so strong to justify 
a one-band model based on localized a'i g orbitals. 

In the following, we test the ability of H e ff with pa- 
rameters fixed as explained above, to explain the ARPES 
experiments. We solve it using a generalization to the 
multiband case of the slave boson treatment of Kotliar 
and Ruckenstein in mean field, which is equivalent to 
the Gutzwiller approximation |23j. The three-band case 
of tig orbitals with onl y o ne relevant spin [2^ and the 
general multiband case |24| were considered before. The 
basic idea is to enlarge the Fock space to include bosonic 
states which correspond to each state in the fermionic 
description. The vacuum state at site i is now repre- 
sented as ej 1 0) , a state with one hole as s\ h\ \0), a 

state with two holes as d\„ , a ,h\ h\ l(T ,\Q), etc., where 



are bosonic operators corresponding 



to empty, singly occupied, doubly occupied sites and so 
on. In our case, with U — * oo, only the first two bosons 
are relevant and from now on we neglect the others in 
the description that follows to simplify it. The boson 
operators should satisfy the constraints 



/ j S iycr S if<? ~ 1) t l i rra .hvy Cr — S^Siyg- (5) 



The first one states that there is only one boson at each 
site and the second that counting states with one hole in 
the fermionic or bosonic representations should lead to 
the same result. In the hopping terms involving different 
sites (i ^ j in Eq.(j2J), the fermion creation operators 
should be accompanied by their respective bosonic part 



l~f<7 



l~f<7 



'1 



eje,- 



" 2 s t e 



1 



c c . 



(6) 



Here the root operators are identical to one if treated 
exactly. They are introduced to reproduce the non- 
interacting limit U = in the mean-field treatment. 

To use the slave boson theory, we first diagonalize 
the on-site part of H\, introducing hole operators h kaa 
with a'i g and e' g symmetry. In the new basis {a} 



{K 



hi i 



e^ 2 )}, the Hamiltonian for holes becomes 



i l c gl/i I c g2 

Hi = + D aa >)h{ aa h haU 

k ca.ot' ,a 

W(fc) = P- 1 {e aal {k) + e' aa ,(k))P 

t a {k) i M (fc) t„,a(fc) 
*a,i(fc) h,i(k) ii, 2 (fc) 

yt a ,2(k) t 2 ,l(k) *2,2(fe) / 

-2 o o^ 

LI,,,' = p- l D aa ,P = D [ 10 

1, 



(7) 



with P being the matrix for the change of basis. In mean- 
field the bosonic operators are replaced by numbers, the 
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values of which are determined self-consistently minimiz- 
ing the total energy under the given constraints. For a 
homogeneous paramagnetic phase we have two different 
slave bosons for singly occupied sites, (si aa ) = s a for a' 
and (s ll 3 a ) = s e for e' . Eqs. © imply 



x=0.3 



x=0.5 



+ 2^ + 4^ = 1; (n) = 2s 2 a + As\ 



(8) 



with (n) the total hole occupancy per site. From Eq. 
the hopping terms become renormalized as i a (fc) — > 
q a t a (k), t e (k) -> q e t e (k) and i oe (fe) -> y/q^qet ae (k), 
where g ra = a;/(l — s„). 

In Fig. |3 we show the resulting dispersion relation for 
electrons (instead of holes to facilitate comparison with 
previous works) for three different values of x near to 
those measured by ARPES, and the corresponding FS. 
A comparison with the case U = in the effective model 
(not shown), for which the band width is near 0.9 eV, 
renders it clear the effect of the band narrowing as a 
consequence of the remaining correlations in the effective 
model. This renormalized bandwidth is stronger as the 
Mott insulating limit x — + is approached. As a conse- 
quence of both, the relatively strong trigonal splitting due 
to the large Co-0 covalency and the band renormaliza- 
tion, the hole pockets near the Brillouin zone boundary 
disappear. We want to stress that if only one of these ef- 
fects were present, these pockets would remain. In agree- 
ment with the most recent ARPES measurements |llj . 
the feature that would give rise to hole-pockets near the 
K points remains under e f with almost no doping depen- 
dence (while the corresponding band bottom shifts from 

100 meV between x = 0.3 and x = 0.5). The FS 

agree with ARPES experiments, showing only a central 
lobe with a hexagonal shape, which is more marked for 
lower dopings. The average amount of a' lg holes on the 
FS is 38% weakly dependent on doping, the a' x charac- 
ter is minimum in the FK direction (~ 31%) and maxi- 
mum in the FM direction (~ 47%). This indicates that 
although only one band crosses the ep, its character is 
strongly mixed and cannot be derived from states of a' lg 
symmetry only. 

In summary, we have derived an effective Hamiltonian 
H e ff to describe the low-energy physics of Na x Co02, 
starting from the ground state of C0O6 clusters, and cal- 
culating the effective hopping between different clusters. 
Previous calculations of the trigonal distortion and fitting 
of the XAS experiment using CoOg clusters fixes the pa- 
rameters of H e ff. This H e ff without any new adjustable 
parameter is able to reproduce the main features of the 
ARPES experiments. 
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FIG. 2: (Color online) Renormalized band-structures for dop- 
ing x = 0.3, 0.5, 0.7 and corresponding FS in red, green, blue 
respectively. Larger FS corresponds to lower doping x. 
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